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Abstract. Gamma-ray bursts (GRBs) have been suggested as possible sources of the high- 
energy neutrino flux recently detected by the IceGube telescope. We revisit the fireball 
emission model and elaborate an analytical prescription to estimate the high-energy neu¬ 
trino prompt emission from pion and kaon decays, assuming that the leading mechanism for 
the neutrino production is lepto-hadronic. To this purpose, we include hadronic, radiative 
and adiabatic cooling effects and discuss their relevance for long- (including high- and low- 
luminosity) and short-duration GRBs. The expected diffuse neutrino background is derived, 
by requiring that the GRB high-energy neutrino counterparts follow up-to-date gamma-ray 
luminosity functions and redshift evolutions of the long and short GRBs. Although dedicated 
stacking searches have been unsuccessful up to now, we find that GRBs could contribute up 
to a few % to the observed IceGube high-energy neutrino flux for sub-PeV energies, assum¬ 
ing that the latter has a diffuse origin. Gamma-ray bursts, especially low-luminosity ones, 
could however be the main sources of the IceGube high-energy neutrino flux in the PeV 
range. While high-luminosity and low-luminosity GRBs have comparable intensities, the 
contribution from the short-duration component is significantly smaller. Our findings con¬ 
firm the most-recent IceGube results on the GRB searches and suggest that larger exposure 
is mandatory to detect high-energy neutrinos from high-luminosity GRBs in the near future. 
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1 Introduction 

Gamma-ray bursts (GRBs) are among the most energetic events in the Universe (see e.g. [1, 2] 
for reviews on the topic) and have been suggested as sources of ultra-high energy cosmic 
rays [3, 4]. In terms of astronomical observations, they are usually divided into two distinct 
groups on the basis of the BATSE bimodal distribution: Long-duration bursts (whose dura¬ 
tion is longer than 2 s) and short-duration bursts (lasting for less than 2 s) [5]. Long duration 
bursts are thought to originate from the collapse of a massive star to a black hole [6], while 
short-duration ones should originate from coalescing neutron stars or black-hole-neutron-star 
mergings [7]. 

It is not excluded that, although with different origin, short- and long-duration bursts 
are driven by the same underlying mechanism; The fireball model [8-10]. According to this 
model, a hot “fireball” of electrons, protons and photons forms. Such fireball is initially 
opaque to radiation, then the hot plasma expands by radiation pressure and particles are 
accelerated making the plasma transparent to radiation and favoring the emission of keV- 
MeV photons [11, 12]. Neutrinos with energies of 0(100) TeV are also expected to be emitted 
from these sources because of lepto-hadronic interactions [8]. Besides the fireball emission 
model [8, 10] on which we will relay in this work, other models have been proposed to explain 
the neutrino production in GRBs, such as the dissipative photosphere model [13] according 
to which the prompt GRB emission occurs near the Thomson scattering photosphere, or 
large-radius magnetic dissipation models, such as the IGMART model that relies on a highly 
magnetised outflow which is dissipated at a radius larger than the internal shock radius [14] . 
See [15] for a common formalism to connect these models and a comparison among the 
expected neutrino fluxes. 
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The IceCube detector, a neutrino telescope made with 5160 optical modules and lo¬ 
cated at the South Pole, could detect neutrinos from GRBs by measuring the Cherenkov 
light from secondary particles produced in the neutrino-nucleon interactions. Over the past 
years, IceCube performed searches for muon neutrinos associated to CRBs, but with negative 
results [16-20]. The IceCube upper limits therefore started to put tight constraints on the 
expected CRB flux and on the theoretical models employed to explain the neutrino emission 
from these sources [21, 22]. The diffuse neutrino emission from CRBs is also widely discussed 
in the literature [23-31] and it has been recently invoked [32-39] as a natural possibility to 
explain the PeV neutrino events discovered from IceCube [40-43] , assuming such events have 
a diffuse origin (see [44] for an overview on the possible astrophysical sources of the IceCube 
PeV neutrinos and references therein). 

In light of the unsuccessful most recent IceCube stacking analysis [43] and of the high- 
energy neutrino flux discovery [40-43], the aim of this work is to provide an up-to-date esti¬ 
mation of the expected high-energy neutrino prompt emission from long- and short-duration 
CRB families within the fireball model [9, 10]. To this purpose, for the first time we present 
an analytically extensive modelling of the neutrino emission from fireballs by including ra¬ 
diative, adiabatic and hadronic cooling processes involving pion and kaon decays. Then, 
we discuss the relevant processes for each CRB family (high-luminosity, low-luminosity and 
short CRBs) and derive an estimation of the diffuse flux by adopting up-to-date gamma-ray 
luminosity functions and evolution of formation rates to fix the normalization of the CRB 
neutrino energy spectrum. 

Our results show that, for average CRB parameters compatible with observations, the 
diffuse neutrino emission from CRBs could contribute up to a few % to the currently observed 
IceCube high-energy neutrino flux in the sub-PeV region while CRBs might be the main 
source of the IceCube flux in the PeV range, assuming that the latter has a diffuse origin. 
Moreover, we find that the neutrino emission within the hreball model is widely compatible 
with the current IceCube limits due to stacked CRBs. 

This manuscript is organised as follows. In Section 2 we introduce constraints coming 
from gamma-ray observations on the luminosity functions and redshift distribution of long- 
duration (divided in low-luminosity and high-luminosity CRBs) and short-duration CRBs. In 
Sec. 3, we analytically model the neutrino emission from py interactions in CRBs and define 
the expected neutrino spectrum from pion and kaon decays including radiative, adiabatic 
and hadronic cooling processes. In Sec. 4, we present our results on the diffuse emission 
from these sources, discussing the astrophysical uncertainties as well as the ones related to 
the CRB model parameters, and compare them with IceCube bounds and future searches. 
Conclusions are presented in Sec. 5. 

2 Observational constraints on gamma-ray bursts 

Camma-ray bursts have been monitored since long time now. They are usually divided into 
two families according to the observation of these transients in photons: Long- and short- 
duration bursts. However, current data are still insufficient to define accurate luminosity 
functions, especially at high redshifts. Sources of errors are, for example, the triggering crite¬ 
ria of the detection instruments that can be responsible for poor estimates in the parameters, 
degeneracies arising from a mixing of the luminosity function and the source rate evolution 
with the redshift, as well as selection effects. Civen such uncertainties, in the following we 
will dehne allowed bands for the distribution of the CRB families in luminosity and redshift 
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Table 1. GRB parameters adopted in the estimation of the diffuse high-energy neutrino flux for 
our canonical model including the astrophysical uncertainties. The local rate (po) is in units of 
Gpc“^ yr“^, the isotropic luminosities (L*, Lmin and Lmax) are expressed in units of 10^^ erg s“^, the 
variability time ty is in s. The best fit parameters a and j3 are employed in the luminosity function 
fits, while a-y and (3^ describe the gamma-ray Band spectrum, and F is the bulk Lorenz factor of the 
jet. 



Po 

T* 

a 

/3 


-^'max 

a-y 

I3y 

T 

ty 

HL-GRB 

0.8 

0.8 

-0.95 

-2.59 

10"^ 

i(F" 

1 

2 

500 

0.1 

HL-GRB 

0.5 

0.48 

-0.13 

-2.42 

10-3 

10^ 

1 

2 

500 

0.1 

LL-GRB 

2000 

5 X 10-4 

-2.3 

-1.27 

1.8 X 10-^^ 

10-3 

1 

2 

5 

100 

LL-GRB 

200 

5 X 10"^ 

-2.3 

-1.27 

1.8 X 10-® 

10-3 

1 

2 

5 

100 

sGRB 

4.6:|- 

2 

-1.94 

-3.0 

0.5 X 10-3 

10 

0.5 

2.25 

650 

0.01 


and characterize the source energetics. Unless otherwise specified, we will distinguish among 
three different reference frames: The GRB reference frame, the jet comoving frame and the 
observer frame. Each physical quantity X will be labelled as X, X', and X in each of these 
frames respectively. 


2.1 Long duration gamma-ray bursts 

Long-duration GRBs are thought to originate from the collapse of a massive star into a 
black hole [6]. They are usually divided in two sub-categories: High-luminosity (HL) and 
low-luminosity (LL) GRBs. 

The source rate evolution as a function of the redshift (z) of the HL component is 
described through a piecewise function [45]: 


Rhl-grb = Po 


(1 + zY for z < Zi, 

(1 + + zY for z > Zi, , 


( 2 . 1 ) 


being pQ the local rate defined as in Table 1, z* = 3.6, a = 2.1, and b = —0.7 [46]. 

The HL-GRB typical luminosities vary in the range [10^®, 10^^] erg s“^. We adopt the 
intrinsic isotropic luminosity function (LF), i.e., corrected for beaming effects and defined in 
the GRB frame, found in [47] on the basis of a selected sample of 175 bursts observed by 
Swift: 


‘f’HL-GRB C)C 


for Liso < L* 
for Liso > T* , 


( 2 . 2 ) 


with best-fit parameters defined in Table 1 [47]. Note that in Table 1 we consider two 
sets for the LF parameters and for the local rate pQ, corresponding to the extremes of the 
accepted band allowed from the data (see [47] for more details). Figure 1 shows the redshift 
distribution of the HL-GRB family in light blue. The uncertainty band is defined by the two 
extreme sets of parameters in Table 1 for HL-GRBs. The distribution of these sources peaks 
at about z ~ z* and it stays roughly constant at higher redshifts. 
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The LL-GRBs are characterised by isotropic equivalent luminosity much smaller than 
that of the HL-GRBs, i.e. erg s“^. The origin of the LL component is not yet clear. 

It is not excluded that the LL-GRBs belong to a distinct population than the HL-GRBs [48] . 
On the other hand, it could be that they have the same origin as HL-GRBs, but are produced 
by failed jets that do not break out of their progenitors. According to the latter hypothesis, 
this could mean that most supernovae generate jets producing LL-GRBs and only few of 
them are able to produce jets powerful enough to break out and produce HL-GRBs [49]. Due 
to their low-luminosity, LL-GRBs have been mainly detected nearby {z < 0.1) and a large 
fraction of their population might be below the detection threshold. 

In order to estimate the LL-GRB total event rate, we assume a supernova-GRB con¬ 
nection, and suppose that the LL-GRB rate should be lower than the one of type Ib/c 
supernovae. Therefore, we define the LL-GRB rate as [23, 50]: 


Rll-grb = Po 


{l + z)P^^ + 


i + zY^^ 


5000 


l + z 


P 3 K- 


1 /k, 


(2.3) 


with K = —10, Pi = 3.4, p 2 = —0.3, ps = —3.5, and we consider two representative values for 
the fraction of SN Ib/c that goes in LL-GRBs [23], i.e., po = 0.01 psN and po = 0.1 psN with 
PSN = 2 X 10^ Gpc“^ yr”^ [51]. The LL-GRB rate is shown as a function of the redshift in 
Fig. 1 (violet band), it peaks at 2 : ~ 1 and it decreases at higher redshifts. Note that it is 
locally much higher than the HL-GRB one. 

The LL-GRB LF is parametrized as in [23, 52]: 


^LL-GRB OC 



(2.4) 


with L*, a and /3 as in Table 1 based on the Swift-HNJ: sample analysed in [52]. Although 
the adopted LF fit describes both the LL and HL-GRBs [52], we use it for the LL-GRB 
component only. 

For both the LL- and HL-GRBs, we assume that the injected (inj) gamma-ray energy 
spectrum is fitted with a Band-spectrum [53]: 


dNY 


OC 



for E'^ < EY , 
for E'^ > EY , 


(2.5) 


with = 1 and Y = 2 as in Table 1. The photon break energy E^^i, is usually expressed as 
a function of the isotropic energy {Eiso) through the so-called “Amati relation” that holds 
for the long-duration GRBs [54, 55]: 


0.1 MeV 


(3.64 ± 0.04) 


Hi, 


7.9 X 10^2 erg 


0.51±0.01 


( 2 . 6 ) 


The isotropic energy of each GRB in the source rest frame can be expressed as a function 
of the isotropic luminosity (Ljso) by combining the “Yonetoku relation” [56] and the Amati 
one [57] for the long-duration GRBs [32]: 


log 


Hi, 


1052 gj.g 


= 1.07 log 


Hi, 


1052 erg/g 


-F (0.66 ±0.54) . 


(2.7) 
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Figure 1. Redshift distribution of the HL-GRBs (light blue), LL-GRBs (violet), and sGRBs (orange). 
Each family is normalised to its local rate po- The LL-GRB local rate is higher than the HL-GRB 
and the sGRB ones, and the sGRB rate quickly decreases for z > 1. 


In the following, we will assume that the Amati and Yonetoku relations hold for both the 
LL- and HL-GRBs, postulating that the two populations share the same emission mechanism 
(see e.g. [48, 58] for a discussion on the topic). 

The photon break energy in the jet frame is related to the same quantity in the GRB 
frame through The normalisation of the photon spectrum, is set by 

assuming that f-y = Eiso/ dEyEy{dNy/dEy). 

2.2 Short-duration gamma-ray bursts 

Short-duration GRBs (sGRB) exhibit typical luminosities similar to those of the HL-GRBs, 
and are believed to originate from neutron-star-neutron-star or neutron-star-black-hole merg¬ 
ers [7]. Therefore, we should expect a delay between the star-formation rate and the merger 
rate due to the spiral-in time [59, 60]. In order to model the sGRB redshift distribution and 
their luminosity function, we follow [61] which selects a sample of non-collapsar GRBs from 
the BATSE, Swift, and Fermi data. 

The sGRB rate is described by the convolution of the ordinary star-formation rate 
with a function /(At) that takes into account the time-delay At of the spiral-in time of the 
binaries [59, 60]. According to [61], the time-delay function that is in better agreement with 
the data is a log-normal distribution that allows to express the sGRB formation with the 
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following effective rate: 


exp[( 2 : — 0.9)/0.39] for z < 
exp[—(z — 0.9)/0.26] for z > z^. 


( 2 . 8 ) 


-RsGRB = 10 Po 


with z* = 0.9. The sGRB rate as a function of the redshift is plotted in Fig. 1 (orange band), 
note as it peaks at about z* and then it decreases quickly. The sGRB rate is locally higher 
than the HL-GRB one. 

The LF is fitted with a broken power law [61]: 


^sGRB oc 


for Liso < T* 
for Liso > T* , 


(2.9) 


with the best ht parameters provided in Table 1. 

Similarly to the long-duration GRBs, sGRBs have a gamma-ray spectrum fitted with 
the Band spectrum (Eq. 3.19). However, we know from observations that the low-energy 
component (i.e., for < E^^i,) is harder for sGRBs than for the long-duration GRBs (see 
the values for in Table 1) and the peak energy is slightly higher [61-63]. 

We assume that relations similar to the Amati and Yonetoku ones hold between 
Eiso, and L^so for the sGRBs. To this purpose, we extrapolate them by fitting the data in 
Fig. 7 of [63] and define the analogous of Eqs. (2.6) and (2.7): 


(lofei) = (iowSiTi) -' 

We suppose that this class of GRBs has shorter variability timescale (t„) than long-duration 
GRBs [64], as reported in Table 1. 


3 Prompt neutrino emission from gamma-ray burst fireballs 


In this Section, we discuss the neutrino production in GRBs through py interactions and 
derive the corresponding neutrino energy distributions. The main reactions that we study 
are: 


p + ^^A^n + 7r~^,p + Tr^ (3-1) 

p -|- 7 —7> K~^ + A/S . 

Pions, kaons and neutrons in turn decay into neutrinos: 


7r+ 








+ 

+ e+ 

Tr~ 








+ Z^e 

1 

+ 

K+ 



+ ^l 

IL 5 

n 


p + e 

+ Z^e • 
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In the following, we will assume that the neutrino contribution from the n decay is negligible 
(see Fig. 2 of [65]) and we will reconstruct the neutrino energy spectrum from the pion and 
kaon decays. 


3.1 Neutrino production from pion decay 

The comoving proton energy threshold to produce a A resonance is E'^ > [(rriAC^)^ — 
(mpC^)^]/(4Fi(,), where ttt-a {'nip) is the A {p) mass. The corresponding neutrino energy 
in the observer frame is [66]: Ey = 0.05 [r/(l + z)]^ {ni\ — mp)/{2E:y), with T the bulk 
Lorenz factor of the jet that we assume constant for sake of simplicity; the numerical factor 
comes from the fact that the average energy fraction transferred from the initial proton to 
the pion is 20% times 1 /4 which comes from the assumption that the 4 leptons in the tt decay 
channel equally share the energy of the pion (i.e., E^, = 0.05 Ep). 

The resultant neutrino spectrum will exhibit a spectral break at: 


Eu,b,tT - 


- 0-77 


1 + z 


. 2\2 


— i'nipC^) 


2E, 


= 1.6 X W'^E. 


6 17-1 


7,b 


7,b,MeV \ 1 _|_ 2 


r 2.5 


GeV 


(3.3) 


where, in the hrst equality, the correction factor of 2 with respect to Ep takes into account 
the fact that the pion production efficiency peaks at higher center of mass energies [65] . The 
pre-factor for neutrinos produced by pion decay is Ott = 0.05 = 0.2 x 1/4 (since 20% 
is the fraction of the proton energy that goes into pion and 1/4 is the fraction of the vr 
energy carried by neutrinos). In order to favor a comparison with the existing literature, 
in the second equality, we express it in terms of typical values for the GRB parameters: 
r2.5 = T/IO^'^ and E-y^b,MeV = ii' 7 ,fe/MeV. 

Above E,y^b,-wj the neutrino spectrum is the same as the proton spectrum, because the 
protons can interact with photons with energies E^^b or smaller, where the number of photons 
is almost constant, i.e., E^dN^/dE^ oc E^f = E!^. Below Ei^^b,n, on the other hand, the 
number of photons with energy E^ that a proton with energy EyjaT^ can interact with 
is suppressed by a factor of {E^/E^^b)^~^i = {Ey^b,-K/, and the resulting neutrino 
spectrum is harder by a factor of Ep'^ ^ with respect to the proton spectrum (that we 
assume to scale as Ep ^). 

If pions produced by the pry interactions decay faster than they cool, the correspondent 
neutrino spectrum is not affected. Otherwise, other neutrino break energies are determined 
by the radiative cooling (rc; synchrotron radiation and inverse Gompton scattering), adiabatic 
cooling (ac), and hadronic cooling (he) processes, being the total cooling timescale defined 
as The neutrino spectrum is modified by an additional factor of 

[1 — exp(—f/mTr/F^/rTr)], where is the pion lifetime and m-,^ is the pion mass. If the cooling 
time scale is much shorter than the decay lifetime in the jet frame (i.e., <C E'^Tj^/ m-,^), 
the additional factor can be approximated as t'^riTr/E'^ t-,^. Below, we shall discuss the three 
relevant cooling processes. 

The hadronic cooling timescale for pions in the jet frame is: 


E'^ 207rc‘^r®mpt^ee 

cann'pAE'^ + z)'^ ’ 


(3.4) 


where c is the speed of light, n'p = Ej{l + z)/iJ'mpiP'V') the jet comoving proton density, 
(T/i = 5 X 10“^® cm“^ [67] the cross section for meson-proton collisions, and AE'^ the energy 
lost by the incident meson in each collision, which is AE'^ = 0.8E'^ [68]. For a jet with kinetic 
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energy Ej, opening angle 9j and variability timescale ty, the volume of an infinitesimal jet 
shell is V' = 2-K6‘jf'jCtjT/{l + z) with fj = 2r^ct^/(l + z) the internal shock radius. We 
then converted the jet energy Ej divided by the duration tj (the jet luminosity Lj = Ej/tj) 
into the isotropic peak luminosity via Lj/(27r6<|) = 0.3Liso/[47r(l + z)^ee], where 0.3 comes 
from the fact that typically 0.3 times the peak luminosity is the luminosity averaged over the 
duration of the burst [32, 69], and eg is the energy fraction carried by the electrons. When 
the hadronic cooling is the dominant cooling process (i.e., t'g ~ Chc)> ^ steepening of the 
neutrino energy spectrum occurs at an energy satisfying = TT^E'^KniTrC^) with Ttt 
the pion lifetime (mass), i.e.. 


E. 


z/, 7 r,hc — 


207rm^mpC^r6tgee ^ ^ ^ ^q13. 

(1 + ^.‘iTT.ahLiso 


2.5^v-2^e 

(1 “1“ '2^)^.£iso,52 


GeV , (3.5) 


where 6^ = 1/4, £'/ = E'^/4, ty -2 = ty/{W~‘^ s), and Ligo ,52 = hso/iW^"^ erg s“^). Since 
is independent of energy, the neutrino spectrum suppression factor due to the hadronic 
cooling is thc^Tr/iKT^r) OC l/E'^. 

The adiabatic cooling time in the jet comoving frame is given by 


t' = 

^7r,ac 


Tc 


(3.6) 


which is energy independent. By solving = Tt^E'^/{ ruT^c^), we find that the neutrino 
spectrum breaks at 


E, 


i/', 7 r,ac 


(1 + zY 


= 2.7 X 10' 


,9^2.5^’'>-2 


{l + zf 


GeV 


(3.7) 


and it steepens with respect to the parent spectrum by 1/iil/. 
The radiative cooling time is given by 


3mic^ 


37rc®mlr®t?, 


4aTm2^;(C/^ + [//) 0.3c7Tm2T;;Zigo(l + 2)2(l + es/ee) 


(3.8) 


where U'^ = B'‘^/{S'k) = 4eBE'-/V' defined in terms of the fraction of the internal energy 
carried by the magnetic field (e^), C// = E'^n'^ = 4£''ee/V',^ mg the electron mass and ut 
the Thomson cross section. Note that oc IjE'^. When t'g ~ and the other cooling 
processes do not play a role at any relevant energies, a break energy in the neutrino spectrum 
is expected for = Tt^E'^/[ rny^c?')-. 


Ey 


b^r 


37rc*m®r®t2 


V 


(l + z)2 [0.3r.,r<7T"ie-^iso(l + es/Ce) 


1/2 


= 1.8x10® 


^2.5^1^,-2 


(l + z)2[Ligo,52(l + eB/ee)]V2 


GeV 


(3.9) 

and the radiative cooling is responsible for a steepening of the correspondent neutrino spec¬ 
trum by a factor oc l/El"^. 


^The numerical factor 4 in the definition of Ug and [7/ comes from the fact that, assnming that the upstream 
material is cold while the shock wave is propagating, the relativistic strong shock transition relations predict 
a post-shock energy density U 2 = IFlime, being r 2 i — 1 the relative Lorenz factor for the internal shock 
model and ni the density before the shock front, and e the energy carried by the particle [2]. 
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If the hadronic cooling occurs at energies lower than the ones for which the radiative 
cooling dominates, the transition between the two cooling mechanisms happens at the neu¬ 
trino break energy corresponding to instead of it is: 




b^T 


Sc^m^cTh 


U,77,VC,2 - 


(1 -I- z) 2QmpmlaT{ee + es) 


= 1.4 X 10'^ 




(1 -|- z)(€e + Cb) 


GeV . 


(3.10) 


Similarly, if the adiabatic cooling occurs before the radiative cooling, the break energy in the 
neutrino spectrum due to the radiative cooling is determined by 


E, 


u,77,rc,3 


_ b7rc^m% =l2x 10^__ 

(1 -k z)2 (1 + eB/ee)UsomlaT (1 + + eB/ee)Liso,b2 


GeV . (3.11) 


In both cases the correspondent neutrino spectrum will be subject to a further !/£'(, steep¬ 
ening with respect to the one affected by hadronic or adiabatic cooling. 

Muons produced by the pion decay (from now on indicated as will in turn originate a 
neutrino energy spectrum, with a first break energy defined similarly to the one in Eq. (3.3), 
but with the pre-factor = 0.05 = 0.2 x 3/4 x 1/3 (where 3/4 is the energy fraction 
transferred from pions to muons, and 1/3 is due to three-body decay of the muon). In terms 
of typical GRB parameters, it can be written as 

= 1.6 X 10‘E-l„,y (it/) 


Other energy breaks in the neutrino spectrum are determined by the radiative cooling 


processes 


hpT 


fj,/ 


■^u,^,vc 


{l + zf 


37rc®m^r®t^ 


11/2 


0.3r/,crTm2Liso(l + es/e 


= 1.3x10^ 


^2.5^11,-2 


(l + z)2[Li,o,52(l + 6B/ee)]V2 

(3.13) 

with bp = 1/3 since E/ = E'p^/3. Note as the break energy Ey^p^c is usually one order of 
magnitude smaller than the pion one (Eq. 3.9), due to difference in the masses and lifetimes 
of the two particles. Similarly to pions, muons are also subject to adiabatic cooling. In this 
case, the correspondent break energies in the neutrino spectrum are 


GeV 




2mu.c^E 

M M— V ^ 3^2 X 10^ 


E, 


i/,^,rc,3 — 


(i + zy 

bpT^ 


^2.5^v,-2 

(1 + ^)2 


GeV , 


bnc^rripty 


(1 + zy (1 -k EB/ee)L\somlaT 


= 5.4 X 10' 


,6 


^15^11,-2 


(1 -k z)^{l -k eB/ee)//iso,52 


(3.14) 

GeX3.15) 


Assuming that E'^ = 4£'i,(l -k z)/T and = 3Eli^(l -k z)/T, the relations among the 
comoving muon and pion lifetimes and the cooling times as from Eqs. (3.4), (3.6), (3.8) 
and the correspondent ones for muons are shown in the top panels of Figs. 2, 3 and 4 for 
typical HL-, LL-GRB and sGRBs at 2 ; = 1. We took Ziso = 10 ^^ erg/s for the HL-GRBs, 
Liso = 10'^® erg/s for the LL-GRB, and Liso = 10^^ erg/s for the sGRBs. For each family, we 
supposed r and ty as in Table 1 and defined Cg = €b = 10“2 for HL-GRBs and sGRBs, while 
we adopted = eb = 10~® for LL-GRBs (see Sec. 3.3). 

For the GRB parameters assumed in Figs. 2, 3 and 4, the hadronic cooling is negligible 
for HL-GRBs and sGRBs. On the other hand, the radiative cooling and the adiabatic one 
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Figure 2. Top panel: Muon and pion lifetimes and cooling times in the jet comoving frame as a 
function of the neutrino energy E^, for a typical HL-GRB with Liso = 10^^ erg/s and z = 1. Bottom 
panel: Muon and kaon lifetimes and cooling times for the same HL-GRB. For the assumed HL-GRB 
parameters, the radiative cooling is always important, while the hadronic cooling is negligible and the 
adiabatic cooling is relevant for muons. 


are always relevant for all three families for both pions and muons. As we will discuss in 
the following, note as such hierarchy among the cooling processes is a function of L^sq and 2 
for other fixed GRB parameters. Therefore it will change within the luminosity and redshift 
range that we will consider for the computation of the diffuse neutrino emission. 

3.2 Neutrino production from kaon decay 

Yet another contribution to the total neutrino spectrum from GRBs originates from kaon 
decays. The resultant neutrino spectrum will have a first break energy coming from the 
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Figure 3. The same as Fig. 2, but for a typical LL-GRB with Ljso = 10“^® erg/s and z = 1. For the 
assumed LL-GRB parameters, the adiabatic and the hadronic cooling are negligible for kaons. 


proton threshold energy for kaon production, similarly to Eq. (3.3): 

(it:;)-21-^- 


(3.16) 


where, for neutrinos directly produced from kaon decay without going through muons, ck = 
0.1 (since 20% is the fraction of the proton energy that goes into K and 1/2 is the fraction 
of K energy carried by neutrinos). In the case of muons originating from kaon decay, from 
now on indicated as /i/c to distinguish them from the ones from muon decay (/Ttt), = 

0.033 = 0.2 X 1/2 X 1/3. In terms of typical parameters, the first break energies for K and 
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Figure 4. The same as Fig. 2, but for a typical sGRB with Liso = 10®^ erg/s and z = 1. For the 
assumed sGRB parameters, the hadronic cooling is negligible. 

fj,K become: 

£^.,b,K = 8.5 X (r^) 

Eu,b,iiK = 2.6 X 10®£'..^GeV . (3.18) 

The resultant neutrino spectrum from kaon decay is affected, similarly to what was 
discussed for pions, from hadronic, adiabatic, and radiative cooling processes. Therefore, we 
should expect breaks in the neutrino spectrum similar to the ones described in Eqs. (3.5), 
(3.7), (3.9), (3.10), and (3.11) according to the cases, but with the substitution m 7 r(r.n-) —)• 
mK{TK)- For example, will be (iiyr/(l+z)^ [47rc®m®r®t^/(0.3r7rCJTmgLiso(l+eB/ee))]G2, 

with dx = 1/2 since F/ = E'^/2. Similarly, the equations in terms of typical parameters 
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can be extrapolated from the ones from pions rescaling m,r('^ 7 r) to mK{TK) and taking into 
account the different multiplicity factors. 

Muons are also produced from kaon decay. The resultant neutrino energy spectrum will 
be determined as the one described for muons from pion decay (see Eqs. 3.13, 3.14 and 3.15). 

The bottom panels of Figs. 2, 3 and 4 show the relations among the muon and kaon 
lifetimes and the cooling times, similarly to the top panels for pions and muons from pion 
decays, assuming that E'j^ = 2Ey{l + z)/T and E'^ = 3E^{1 + z)/T. The hadronic cooling is 
negligible, while the adiabatic cooling is relevant for muons for the model parameters plotted 
in these figures. Note that the break energies in the neutrino spectrum due to the kaon 
radiative cooling occur at higher energies than the ones due to pion radiative cooling given 
the differences in the rest-mass and lifetimes of the two parent particles [70, 71]. We neglect 
here the contribution to the neutrino flux coming from , see Ref. [38] for a dedicated 
discussion. 


3.3 Neutrino energy spectra 

As discussed in the previous section, the neutrino energy spectrum for one flavor [u + v) 
resulting from vr, AT, and /r decays will be a broken power law derived from the parent proton 
spectrum (that we assume is proportional to with further breaks defined according 

to the hierarchy of the cooling processes. For example, when the pion hadronic cooling is 
negligible, the adiabatic cooling is relevant for muons only, and E'^^ 3 < 

the resultant neutrino energy spectrum produced from muons from pion decay will 
be: 


dN^ 
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(3.19) 


For each decay channel i, the neutrino energy spectrum is normalized in terms of the 
total photon fluence by generalizing the expression proposed in Refs. [10, 72]: 
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dNy \ 
dEy ) 


(3.20) 


In our numerical computations within the canonical model, we assume the gamma-ray- 
proton luminosity ratio hyp = Liso/Lp = 10“^ for HL-GRBs and sGRBs as suggested by joint 
analysis of high-energy neutrino and ultrahigh-energy cosmic ray (UHECR) data [28, 30] and 
under the assumption that sGRBs behave similarly to HL-GRBs; while we adopt hyp = 10“^ 
for the LL-GRB family assuming that the 7 production is suppressed in these GRBs with 
respect to the other two GRB families for the same Lp^. We will discuss in detail the 

^Note as, assuming Lj = Lp + Le +Lb and ^ tB, one has that h^p ~ ee/(l — Ce — es) ~ £e. This justihes 
the choice of the numerical values of Ce and es introduced in Sec. 3 and adopted trough the whole paper. 
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dependence of the diffuse background from this parameter in Sec. 4.2. For the case of 
neutrinos produced from the pion decay, = 0.97/8 (the coefficient 1/8 = 1/4 x 1/2 since 
each neutrino takes about 1 /4 of the pion energy and 1 /2 of the produced pions are charged 
pions, while 97% is the probability that pions are produced from one py interaction [73, 74]) 
and ~ ,r,second/-E'iy, 7 r,first)/ln(F;p,max/-E'p,min) with first the minimum neutrino 

break energy and second the second neutrino break energy, both defined according to the 
hierarchy of the break energies determined by the cooling processes and Note that, 

for each GRB family and fixed parameters, such energy hierarchy varies as a function of Liso 
and z. The minimum proton energy is i?p,min = rmpC^/(l + z) while the maximum proton 
energy is £'p,max and it is defined as the energy when the acceleration time t/^acc = ^p/{B'ec) 
equals to t/ ,, = min(t/being t/^gync = (SnipC^STr )/the synchrotron 
cooling time and tp^yn = r'Jc the dynamical time scale. For p produced from vr decay, 

— 0.97/3 X 3/8 — A/r and hp^p^ — ln(F/j/p gecond/T'jy^p^first)/hi(F/pjmax/T'p^min) defined 
analogously to hp^p. 

In the case of neutrino production from kaon decay, Nk = 0.03 x 0.63/2 where the nu¬ 
merical coefficient 1 /2 comes from the fact that each neutrino takes about 1 /2 of the kaon’s en¬ 
ergy, 0.63 is the probability that K decays in neutrinos, and 3% the probability that kaons are 
produced from one py interaction [73, 74]. The term hp^K ^ ln(£'j.,if,second/-E'iy,iC,first)/ ln(F;p,max/T'p,min) 
is defined similarly as hp^T^. For muons generated from kaon decay, we have Np^^ = 0.03 x 
0.63/6. 

The second term in Eq. (3.20), / = [1 — (1 — (Xp))^^"^], represents the fraction of the 
proton energy that goes to pion production with {xp) — 0.2 the average fraction of energy 
transferred from protons to pions or kaons per py interaction [10, 15]. The py optical depth, 

Tpry, is defined following [ 8 , 9, 15]: 


Tpn — 


Tt 


= 0.6 




■pi 


1052 erg/g Y102-5 


-2 


lO^^cm 


(3.21) 


with (.p^ the mean free path for py interactions. Note that with respect to the numerical 
fireball model revised in [10], we do include their correction term Cs (originally introduced 
in [72]), but do not implement the factor due to the pion efficiency (see [10] for details). 
The latter is responsible for a further reduction of the expected neutrino flux with respect 
to the one obtained from our Eq. (3.20), as shown in the left panel of Fig. 1 of [10]. We 
also neglect multi-pion processes that should be instead responsible for an enhancement of 
the neutrino flux (see right panel of Fig. 1 of Ref. [10], comparison between their RFC and 
NFC models). Therefore our analytical model is still an approximation with respect to the 
full numerical treatment, but it gives us results in good agreement with the revised fireball 
model of [10] for the same choice of the initial parameters. 

In order to consistently normalize the neutrino spectrum as in Eq. (3.20), we need to 
select the CRB parameters in such a way to have optically thin media [8]. Therefore, we 
consider a test photon with energy = 100 MeV and define the yy optical depth 


_ C' _ (1 + ^)^o'T0-3TisoE-.y^t 
Ti^-y 1287rr®c^t„(meC^)^ ’ 

with the mean free path for pair production i~y ~ (aTUyEy^.)/[16{meC^)‘^] and Lave = 
0.3Liso = ‘i'jrfjT'^cUy [8, 9]. In our canonical CRB model, we fix to the average value 
preferred by observations [23, 64, 75] and define F as the minimum Lorenz factor required 
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to avoid high pair-production optical depth that would inhibit the gamma-ray emission (i.e., 
7'77 ^1) for the whole [L[so, z] parameter space of the HL-GRB and sGRB families considered 
in Sec. 4.1. The Lorenz factors T defined in this way are reported in Table 1. Note that 
belongs to the upper tail of the observed photon energy spectrum and therefore it allows 
us to define an upper limit for the optical depth for pair production for the assigned GRB 
parameters. The < 1 condition does not apply to LL-GRBs as the average photon energy 
is lower than the pair production threshold (see Sec. 2). Variations of the expected diffuse 
neutrino background as a function of ty and T will be discussed in Sec. 4.2. 

The observed neutrino spectrum from a single source at redshift z is dehned as 


F,{E,) 


(l + z)3 UNA 

inTdliz) [dE'J ’ 


(3.23) 


with E' = E{1 + z)/T and diiz) the luminosity distance [76] computed assuming a flat 
ACDM cosmology with 0,^ = 0.32, = 0.68 and Hq = 67 km s“^ Mpc“^ for the Hubble 

constant [77]. The total Ue and neutrino spectra from pion decay at the source and without 
flavor oscillation are: ((iiVi.^/dR'(,)inj, 7 r = {dNy/dE'^)^^ and {dNy^/dE^A^T, = {dNy/dEl)^^ + 
{dNy/dE'^)-,^ and similarly for kaons. Note that no r neutrinos are produced. 

Neglecting flavor oscillations, Fig. 5 (top panel) shows the neutrino energy spectra 
coming from pion, kaon and muon decays as a function of the energy for a typical HL-GRB 
and normalized as in Eq. (3.23) for a source at z = 1. The neutrino spectrum coming from 
muon decay from pions dashed blue line) exhibits four breaks, the first one corresponds 
to the first energy break (Ej^,b,^), the second is due to the radiative cooling {Ey^^Ai the third 
to the adiabatic cooling of muons {E^^^A the fourth occurs at the same energy of the 
second break of the pion spectrum (blue line) and it is indeed due to the radiative cooling of 
the parent pion see Fig. 2 for comparison). For kaons (magenta lines) only the first 

break energy is relevant, while the neutrino spectrum coming from muon decay from kaons 
(/XX, magenta dashed line) is similar to the /Xtt one. The bottom panel of Fig. 5 shows the 
resultant Ve (in black) and (in red) neutrino energy spectra from pion and kaon decays 
and as a function of the energy without flavor oscillations for a HL-GRB at z = 1. Note 
as the cutoff in the neutrino energy spectrum corresponding to the maximum proton energy 
appears for energies slightly larger than the ones studied here. 


3.4 Neutrino flavor oscillations 

While neutrinos travel to reach the Earth they are subject to flavor oscillations and therefore 
the observed flux will be (see e.g., Appendix B of [44]): 


dNy \ 

dE'y) 


^sin2(200) 



(3.24) 


with sin^(20Q) ~ 8/9. Note that while iZy are not produced at the source, the three neutrino 
flavors are equally abundant after flavor oscillations. 

Figure 6 shows the expected E^E^c (computed through Eqs. 3.23 and 3.24) as a function 
of the energy for a typical HL-GRB (Liso = 10®^ erg s“^, blue line), LL-GRB (Liso = 
1048 gj.g g-i^ green line), and for a sGRB (Liso = 10^^ erg s“^, red line) located at z = 1 and 
including contributions from pion, kaon and muon decays. Note as the HL-GRBs give the 
highest flux as expected, while the LL-GRB and the sGRB neutrino spectra are two-three 
orders of magnitude smaller than the HL-GRB one for the adopted input parameters. Note 
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Figure 5. Predicted E‘^Fv{E^) for a typical HL-GRB {Liso = 10®^ erg s“^, 2 = 1) without flavor 
oscillations. Top: Neutrino fluence from tt (blue line) and K (magenta line) decays as well as from 
^ from pion decay (/r^, dashed blue line) and ^ from kaon decay dashed magenta line). Bot¬ 
tom: ElF^^{Ey) (black line) and E^F^^{dE^) (red line) for a typical HL-GRB and without flavor 
oscillations. 

as by adopting the analytical prescription developed in Sec. 3, we find that our estimation 
of the neutrino flux from GRBs gave results close to the ones obtained adopting numerical 
routines in [10, 20, 74], by adopting their same GRB inputs. 
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Figure 6. Predicted E1F''i^{E'^) for a typical HL-GRB (Ljso = 10®^ erg s“^), LL-GRB (Ljso = 
1048 gj-g s“^), and sGRB [Liso = 10®^ erg s“^) at z = 1 with flavor oscillations included. The BL- 
GRBs exhibit the highest flux and the kaon contribution affects the high-energy tail of the spectra in 
all cases. 


4 High-energy diffuse neutrino background from gamma-ray bursts 


In this section, we present our results on the high-energy diffuse neutrino background from 
GRB fireballs. We first discuss the expected neutrino background within the canonical model 
in terms of the astrophysical uncertainties on the local GRB rates and luminosity functions 
(see Table 1), then we study the dependence of the high-energy diffuse neutrino flux from 
the model parameters for each GRB family (see Table 2). 


4.1 Expected diffuse background and uncertainties on the local rate and lumi¬ 
nosity function of each GRB family 

The diffuse neutrino intensity from each GRB component (X) can be defined in terms of 
the gamma-ray luminosity function, through $x(Tiso)dTiso = ^x{Ly)dLi, with the LF 
introduced in Sec. 2 (normalized to unity after integration over luminosity): 


^^max pLxnayi 

IxiE,) = dz dU 

•^^min *^-^min 


z)^ EQa 


Rx{z)^x{Lii 


dNu, 

dE' 


(4.1) 


In the numerical computation of the neutrino background, we assume Zmin = 0 and Zmax = 
11, Liso G [Tmin,Tmax] with I/min and Tmax defined as in Table 1 for each family X, and 
E'^ = Eu{l -b z)/T. Note as the chosen values for ty and T (Table 1) should guarantee us to 
extrapolate an average description of the whole GRB population. However, our estimation 
of the diffuse neutrino emission from GRBs also depends on parameters such as Ce, cb, T and 
h^p that are currently poorly constrained from observations (see discussion in Sec. 4.2) and 
should therefore be considered with caution. 
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Figure 7. Diffuse intensity as a function of the neutrino energy after flavor oscillations for the 
HL-GRB (blue band), LL-GRB (violet band) and sGRB (orange band) families. The bands represent 
uncertainties related to the luminosity functions and local rates (Table 1), whereas all the other 
GRB parameters are fixed to the canonical values. The best fit estimation of the high-energy diffuse 
neutrino flux as in [43] is plotted in light blue, while the blue dot (IC-GRB) marks the upper limit of 
the GRB diffuse neutrino flux from the IceCube Collaboration [20] . The diffuse neutrino background 
from GRB fireballs is smaller than the observed high-energy IceCube neutrino flux in the sub-PeV 
energy range and it scales differently as a function of the neutrino energy. 


For each population X, we implement the analytical recipe described in Sec. 3 and auto¬ 
matically define the neutrino energy spectrum according to the specific hierarchy among the 
different cooling processes for each {Liso^z). Note as for luminosities and redshifts different 
than the ones adopted in Figs. 2, 3 and 4, the hierarchy among the cooling times changes. 
For example, we find that the adiabatic cooling becomes relevant for pions and kaons when 
Liso is on the lower tail of the studied luminosity interval for all the three GRB families. 

We do not include HL-GRBs and sGRBs whose parameters (Liso, z) violate the condition 
"^77 < 1 (Eq. 3.22) in our calculations. However, for the assumed input parameters, > 1 
is realized only for sources with z > 7 and with luminosities at the upper extreme of their 
interval. Therefore, our computation might underestimate the expected diffuse flux only by 
a few % since the diffuse neutrino flux is not affected from sources at z > 7. 

Figure 7 shows the diffuse high-energy neutrino intensity for the HL-GRB (light-blue 
band), LL-GRB (violet band) and sGRB (orange band) components as a function of the 
neutrino energy. Each band takes into account the uncertainty due to the LF determination 
as from Table 1. 
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Our results should be compared with the recent IceCube discovery of high-energy neutri¬ 
nos [43] (blue line in Fig. 7), whose sources are still unknown as well as with the unsuccessful! 
searches on GRBs from the IceCube telescope [16-20]. The total estimated diffuse flux from 
the GRB prompt emission can be as large as 2 x 10“® GeV cm“^ s“^ sr“^ at 10® GeV and it 
is therefore slightly lower than the best ht of the high-energy neutrino flux detected with Ice- 
Cube (blue line in Fig. 7) between 25 TeV and 1.4 PeV. Our estimated total diffuse emission 
is smaller than the IceCube flux at lower energies and it scales differently as a function of the 
neutrino energy. This implies that, assuming that the hreball model properly describes the 
GRB neutrino emission, GRBs cannot be the major contributors to the observed IceCube 
flux for the sub-PeV region. 

IceCube results based on the monitoring of 505 observed GRBs are presented in [20] 
with four years of data. Extrapolating from the high-energy neutrino flux recently detected, 
it is estimated in [20] that the GRB diffuse emission at 100 TeV should be smaller than 
2 X 10“^® GeV cm“^ s“^ sr“^ (blue dot in Fig. 7, IC-GRB). Although it might be that such 
extrapolation on the GRB emission underestimates the real GRB diffuse flux, as IceCube is 
sensitive to the first break energy that occurs for E^, > 100 TeV (see Fig. 6), such upper 
bound is very close to our results for the HL-GRB family. In fact our estimated high-energy 
neutrino intensity at Ei, = 100 TeV is 8x 10“^^ GeV cm“^ s“^ sr“^ for HL-GRBs. Comparing 
the contribution from the HL-GRB component with the exclusion limits presented in Fig. 1 
of [20], we conclude that our results are compatible with current observations and in rough 
agreement with the updated computation of the Waxman-Bahcall flux [8] proposed by the 
IceCube Collaboration (see Fig. 1 of Ref. [20]). Rescaling the HL-GRB flux in Fig. 6 to 
the proton-to-photon luminosity ratio adopted by the IceCube Collaboration (hp^ = 0.1), at 
5 X 10® GeV, our predicted stacking flux for 505 sources is 9 x 10“^ GeV cm“^, which is 
roughly in agreement with the prediction presented in Fig. 2 of [20] for the standard GRB 
model. ^ 

The total expected neutrino emission from GRBs might be further enhanced by GRBs 
that do not trigger the detector because of their low flux or because they are dark in gamma- 
rays [23, 70, 78]. However, such contribution is not discussed in this work. We also estimated 
that the contribution from pp interactions should be subleading with respect to the one from 
p'j interactions in the energy range of relevance. We expect that the overall normalization 
of the expected neutrino background from GRBs should be independent from the employed 
GRB model for fixed GRB parameters (see, e.g. Fig. 1 of [15]), although peaking in slightly 
different energy intervals. 

4.2 Expected diffuse background and uncertainties on the jet parameters 

Up to now, we relied on the canonical models for each GRB family. In this Section, we 
adopt the GRB luminosity functions and local rates corresponding to the upper limits of the 
expected diffuse backgrounds plotted in Fig. 7 and discuss how our estimation depends on 
the photon-to-proton luminosity ratio h^p, the bulk Lorenz factor F and the variability time 


- Dependence on the photon-to-proton luminosity ratio h-^p. In our canonical model we 
assumed h-^p = 10“^ for the HL-GRB and the sGRB components. Such a choice is 

®Note that we focus on the HL-GRB family in the comparison with the IceCube results on stacking searches 
as those are the sources for which there is higher probability to have a photon-neutrino correlation adopted 
in the IceCube procedure for the GRB discrimination. 
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Table 2. Variability range of the GRB model parameters adopted in the estimation of the diffuse 
high-energy neutrino flux. The variability time is expressed in s. 



r . 
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consistent with a coherent picture of neutrinos, gamma-rays and cosmic rays assuming 
that GRBs are main sources of the observed UHECR flux (see, e.g., [10, 28-30] for 
dedicated discussions). Current data suggest lower bounds for this parameter, but 
its precise value has not yet been fixed. Variations of h^p correspond to an energy- 
independent scaling of the neutrino flux shown in Fig. 7 as from Eq. (3.20). 

Assuming that the cosmic ray energy budget of the HL-GRB and LL-GRB populations 
are comparable, h^p for the LL-GRB component should be roughly three orders of 
magnitude smaller than the HL-GRB one (because of the difference in Liso of the two 
populations). However, h^p = 10“^ would give a diffuse neutrino flux above the current 
IceCube observed flux. We therefore chose our canonical value in such a way to boost 
the neutrino emission with respect to the photon one without violating the IceCube 
current bounds [43]. We stress, however, that this parameter is currently unconstrained 
and scalings of the LL-GRB intensity with respect to the one presented here are not 
excluded yet. See also discussions in Refs. [25, 34, 39, 75]. 

- Dependence on the bulk Lorenz factor L. As mentioned in Sec. 3.3, we fixed ty of 
each GRB family to the observed values in the canonical model and determined L to 
guarantee optically thin GRBs for the whole (Liso, z) parameter space. Here, we loosen 
the < 1 constraint and study how the diffuse emission of each GRB family varies 
as a function of L. For the HL-GRB and sGRB components, we define Lmin and Lmax 
inspired from Fermi data [75, 79, 80], as reported in Table 2. The T factor of LL-GRBs 
is poorly constrained due to the scarce statistics collected on these sources up to now. 
However, T of a few seems to be favored (see, e.g., [75, 81] and references therein). We 
therefore consider Tmin = 2 and Tmax = 20. The adopted intervals in T for each GRB 
family are summarised in Table 2. 

Figure 8 (top panel) shows the correspondent diffuse GRB intensity as a function of the 
T parameter. In general, the expected diffuse neutrino intensity increases as T decreases. 
Besides the intensity normalization, T also affects the neutrino break energies as from 
Sec. 3. Note that, by adopting a wide range of variability for the HL-GRB component, 
the diffuse intensity could even be comparable with the IceCube high-energy neutrino 
flux. However, an average T = 100 for the HL-GRB population does not guarantee 
optically thin sources for any L\so and besides being disfavored from the most recent 
IceCube results [20]. 

- Dependence on the variability time ty. The variability time ty has been fixed in the 

canonical model to the average value preferred from observations [23, 64, 75]. Here, 
we study how the diffuse intensity changes for a minimum variability time and 
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a maximum one (t,;,max) defined in Table 2 as suggested from recent data [23, 64, 75]. 
Note that we consider one order of magnitude smaller than the one defined for 

example in Fig. 4 of Ref. [64] assuming that current experimental sensitivity might not 
allow to detect even smaller variability times. 

Figure 8 (bottom panel) shows the correspondent diffuse GRB intensity as a function 
of the variability time ty. Similarly to F, the variability time affects the normalization 
and the neutrino break energies as from Sec. 3. 

In all the studied scenarios, GRBs cannot explain the total detected high-energy Ice- 
Gube neutrino flux in the sub-PeV energy range. However, for certain choices of the model 
parameters, the diffuse neutrino emission from GRBs can reach the IceGube band around 
PeV energies. 

5 Discussion and conclusions 

Gamma-ray bursts (GRBs) are considered to be high-energy neutrino emitters and candidate 
sources of the ultra-high energy cosmic rays. We elaborate an analytical model to compute 
the neutrino flux from GRBs based on the fireball picture, by including the pion, kaon and 
muon decay contributions and taking into account radiative, adiabatic and hadronic cooling 
processes. 

We revisited the hreball model for the GRB neutrino emission [8, 10], providing a simple 
analytical recipe to compute the GRB neutrino flux. For the first time, we included the kaon 
and muon contribution other than the pion one and discussed the role of adiabatic, hadronic 
and radiative coolings. Moreover, for the first time, we investigated the relevant cooling pro¬ 
cesses for different GRB families, i.e., long- (divided in low-luminosity and high-luminosity) 
and short-duration bursts and conclude that while the hadronic cooling is negligible for the 
HL-GRB and sGRB families, the radiative cooling is always relevant, and the adiabatic one 
is not negligible for muons, and for pions of LL-GRBs and sGRBs. 

By adopting up-to-date luminosity functions corrected for beaming effects, we extrapo¬ 
lated the expected diffuse neutrino flux from the long- (divided in low- and high-luminosity) 
and short duration bursts. Assuming that each burst has typical parameters inferred from 
observations and luminosity varying with the redshift as prescribed from the luminosity func¬ 
tions in Sec. 2, we found that the estimated diffuse background intensity could be as large 
as 2 X 10“® GeV cm“^ s“^ sr“^. The low-luminosity GRBs, being the most abundant ones, 
seem to dominate the overall diffuse intensity in the sub-PeV region, while high-luminosity 
GRBs are the major source at larger energies. Such conclusions are also supported from 
variations of the model parameters within the range allowed from observations. We conclude 
that GRBs do not appear to be the leading sources originating the observed IceGube neutrino 
flux in the sub-PeV region, if the latter is interpreted in terms of unresolved sources. Very 
recently, the possibility that the PeV neutrinos originate by GRB cosmic rays interacting 
with the interstellar gas while propagating in the host galaxy has also been investigated [35] . 
However, also such an estimation appears to predict an insufficient neutrino flux. 

Previous estimates of the diffuse neutrino emission from the high-luminosity GRBs [22, 
24, 33] found a total diffuse flux from GRBs of the order of 10“® GeV cm“^ s“^ sr“^ depending 
on the luminosity function and redshift evolution adopted. Our results are in agreement with 
the more conservative ones presented in Ref. [32], where an analysis on triggered and un¬ 
triggered sources has been conducted, normalizing the triggered GRBs through the ones 
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Figure 8. Diffuse intensity as a function of the neutrino energy after flavor oscillations for the 
HL-GRB (blue), LL-GRB (green) and sGRB (red) families, for different values of F (top panel) and 
ty (bottom panel) as from Table 2. The best fit estimation of the high-energy diffuse neutrino flux as 
in [43] is plotted in light blue. 


detected with Fermi/GJiM. In agreement with Refs. [24, 25], we conclude that the neutrino 
flux from the low-luminosity GRBs could be comparable to the one from the high-luminosity 
GRBs and such GRBs could be the main contributors to the IceGube high-energy neutrino 
flux around E^, ~ PeV. 

Our results, although should be considered with caution given the uncertainties on the 
GRB models for the neutrino emission and their parameters, suggest that larger exposure 
is required to discriminate neutrinos from high-luminosity GRBs in forthcoming stacking 
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searches. Moreover, we find that the GRBs do not appear to be the main sources of the 
sub-PeV neutrino flux observed by IceCube. However, if the high-energy diffuse neutrino 
flux resuits from the superposition of severai unresoived sources, for exampie comparing the 
contribution from starburst gaiaxies (see Fig. 4 of [82]) with the one from high-iuminosity and 
iow iuminosity GRBs, it appears that whiie starbursts couid expiain the iow-energy taii of 
the IceGube flux up to 0.5 PeV, the GRBs couid be responsibie for the high-energy taii of the 
neutrino spectrum, if this has a cutoff at energies iarger than a few PeV. Such a hypothesis 
couid be eventuaiiy tested in the coming years in iight of the increasing IceCube statistics. 
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